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Abstract 

Convection structures in binary fluid mixtures are investigated for positive 
Soret coupling in the driving regime where solutal and thermal contributions 
to the buoyancy forces compete. Bifurcation properties of stable and unsta- 
ble stationary square, roll, and crossroll (CR) structures and the oscillatory 
competition between rolls and squares are determined numerically as a func- 
tion of fluid parameters. A novel type of subharmonic bifurcation cascade 
(SC) where the oscillation period grows in integer steps as n^j is found and 
elucidated to be an entrainment process. 

PACS: 47.20.Bp, 47.20.Lz, 47.54.+r 
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Many nonlinear dissipative systems that are driven away from equilibrium by a station- 
ary external forcing respond by spontaneous macroscopic selforganization in time periodic 
oscillating patterns, e.g., in the form of propagating waves but also in standing structures 
0. A system that, depending on parameters, shows both forms is convection in binary fluid 
mixtures confined between horizontal plates that impose a vertical temperature gradient: 
This forcing generates for negative thermodiffusive Soret coupling ip between temperature 
and concentration fields propagating convection waves On the other hand, for ip > 
periodic pulsations of a standing square pattern have been seen such that the intensity 
oscillations are global and homogeneous along the lines of the squares. 

Here we report that and we elucidate how and why such pattern oscillations undergo a 
SC in which the period r does not double but increases as n^j with n = 1,2,3,4, .... This 
peculiar dynamical behavior reflects a novel bifurcation pathway of global pattern oscillations 
that has not been reported before. To explain it we have to keep track of the bifurcation 
properties of three different stationary convective structures (squares, rolls, and CRs as 
explained below). These three states "organize" the dynamical behavior of the system, i.e., 
of the oscillatory convection. They allow to understand the SC as an entrainment process 
and to identify the integer n in the SC as the number of windings around a stable CR state 
of increasing attraction which at the end of the SC becomes sufficiently strong to quench 
the oscillations. 

We have solved the full 3D Oberbeck-Boussinesq hydrodynamic field equations for con- 
vection in binary fluids with two different numerical methods: a finite-difference code 
augmented by a fast multigrid pressure iteration scheme 0] and a Galerkin expansion which 
also yields unstable 3D solutions. Top and bottom boundaries were rigid, impermeable, 
perfectly heat conducting plates. Using a 2 x 2 integration domain in x, y with periodic 
boundary conditions the wavelength A = 2 || was fixed to the experimentally observed 
one @-44|. A few simulations were done for 16 x 16 systems with rigid, impermeable, heat 
insulating side walls to check the robustness of the patterns. We now discuss our findings. 

When the thermal stress, i. e., the Rayleigh number is increased the primary stable 
convective structure for positive Soret coupling, ip > 0, has the form of stationary squares 
for a wide range of Lewis numbers L and Prandtl numbers a as predicted by bifurcation 
theory ||. Then, convection in the form of pure rolls is an unstable solution. Close to the 
convective threshold the square structures can roughly be viewed to consist of a superposition 
of two orthogonal sets of straight parallel stationary convective rolls. We shall call them 
rc-rolls with mode or flow intensity X and wavevector oriented in x-direction and y-rolls, 



respectively, with mode intensity Y and wavevector in ^-direction |10|. In square convection 



the intensities of the two constitutive rolls are the same, X$ = Y$. For the unstably 
coexisting pure a;- roll solution one has Xr > and Yr = 0. The symmetry degenerate pure 
y-roll solution has Xr = and Yr > 0. 

We found that the stationary solution branches for squares and rolls are connected by 
branches of stationary CRs (lines with triangles in Fig. p]) that bifurcate out of the square 
branch and merge with increasing r with the roll bifurcation branch. Bifurcation properties 
of these states were unknown so far despite experimental hints for their existence fl2|,|TT 



In the CR state the two sets of x-rolls and ?/-rolls superimpose with different intensities, 
e. g., Xqr > Ycr > but equal wave number k x = k y . These structures differ from the 
rectangular {k x ^ k y ) CR patterns that grow in pure fluid convection at large r out of 
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rolls ||12|| . The transition from CRs to pure rolls is marked by the vanishing of Ycr- In 
Fig. |l| we identify X C r by upwards (Y C r by downwards) pointing triangles. The CR state 
with interchanged intensities and dominating y-rolls is symmetry degenerate with the former 
and corresponds to a convective structure that is rotated by 90°. In addition to the obvious 
magnitude relations of mode intensities in Fig. ^ like, e.g., < Ycr <Y$ = X$ < Xcr < Xr 
we found for each r shown there that 2Xs ^ Xcr + Ycr ^ Xr, Ns ^ Ncr ^ Nr, and 
wr < wcr < w s ~ V%wr [1131. Here the Nusselt number N behaves as N — 1 oc (X + Y) 
while the vertical flow amplitude w oc (\fX + \/Y). In the above relations the CR quantities 
enter only where they exist, e.g., Ncr starts out of the square state with N C r = N s , grows 
with increasing r, and ends with Ncr = Nr in the roll state. 

For not too large Lewis numbers, i.e., for liquid mixtures the competition between square 
and pure roll convection leads to the appearance of pulsations of X(t) and Y(t) |2|-fl- These 
oscillations grow in a supercritical Hopf bifurcation out of the square state thereby rendering 
the latter unstable. Fig. |2| shows for several Rayleigh numbers that X{t) and Y(t) oscillate 
in opposite phase around a common mean value given by the unstable square state (□ in 
Fig. 0) such that always X(t) +Y(t) = 2X$- The a;-roll intensity X(t) of the pulsating 
pattern grows or decreases on cost of the y-Toll intensity Y(t), however, without ever going 
to zero. Thus the two roll sets never die out completely or reverse their turning direction 
during the oscillations. 

In Fig. ^ we show grey-scale snapshots of the topview shadowgraph intensity distri- 
bution I(x,y) taken at times marked in Fig. |^ by A (X minimal; y-rolls dominant), B 
(X = Y; squares), and C (X maximal; x-rolls dominant). After half an oscillating pe- 
riod the pulsating convective structure appears to be rotated by 90° in the x — y plane 
since T[x,y,z\t + ~) = J-{y y —x,z;t) holds for the fields T = 5T and 5C of temperature 
and concentration, respectively. Furthermore they show always the mirror glide symmetry 
T{x, y, z; t) — — J 7 (x+^, 2/+f, 1 — z; t) Qj. The intensity distribution I(x, y) was calculated 
with the formula |I9 

J = I T + I c ~ r {pi + d 2 y ) £ dz (5T + b5C). (1) 

The weighting factor b = - 0.43 is for the parameter combination ip = 0.23, L = 0.0045, a = 
27 that can be realized by ethanol- water mixtures ||. The shadowgraph intensity profile 
along the dotted line in Fig. ^ B is shown in the bottom right corner together with the contri- 
butions I T and 1° to it from the temperature field and the concentration field, respectively. 
Narrow concentration plumes caused by the smallness of L generate a narrow central bright 
stripe in B at the upflow location x = 0.5 and a peak in / on top of the smooth contribution 
from ST. This characteristic structure should be easily accessible in experiments. 

Close to the Hopf bifurcation the oscillations are harmonic and of small amplitude 
(Fig. §a). With increasing thermal driving r the frequency decreases roughly linearly 
(for example, in Fig. [l| the Hopf frequency is ujr = 0.251 and u(r = 1.173) = 0.0586). Fur- 
thermore, and more importantly, with increasing amplitude the oscillation becomes more 
anharmonic and relaxational in Fig. 0b, c; see also Ref. 0|§: While the system rapidly 
sweeps through the square state (□) it spends more and more time in the vicinity of the 
roll state (o). The change from harmonic in Fig. to strongly relaxational oscillations in 
Fig. 0c is also documented in the right column of Fig. 0. There we show phase space plots 
X, Y versus X, Y associated with the time histories of X(t), Y(t) in the left column. 
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At larger r the system gets attracted into one of the CR fixed points that have become 
stable shortly below the r-interval marked SC in Fig. [l]. In this interval we have observed a 
novel subharmonic bifurcation cascade in which the CR attractors entrain the oscillations: 
First the CR attractors deform the phase trajectory (Fig. ||]c). Then with increasing r the 
oscillations execute an increasing number of windings around the CR states (triangles). In 
Figs. |2]c - f the winding number around a CR fixed point increases from n = 1 to n = 4 and 
the period r of the oscillations increases from — to in integer steps. This increase of the 
winding number in the SC continues beyond n = 4; we have found also n = 5. However, the 
control parameter interval 5r n for an n-cycle becomes so narrow — 5r n = (1.6, 0.6, 0.2) 
for n = (2, 3, 4) — that our numerical resources were not sufficient to resolve the SC further. 
But we think that the SC is a robust, experimentally accessible phenomenon; Ref. [[| contains 
a hint for a 2-cycle. Beyond the SC interval in Fig. [1] the system gets attracted into one of 
the CR fixed points. The transition between the oscillations and the CR state is slightly 
hysteretic. Upon reducing r the system remains in the CR state until below the SC interval 
and then a transition to an n = 1 oscillation occurs. So in a small r- interval there is bistable 
coexistence of oscillations and stationary CRs. 

In Fig. |] we present phase diagrams in 2D planes of the r — ip — L — a parameter space. 
Consider first the L — r plane of Fig. |]a. Upon crossing the bifurcation threshold of the 
conductive state squares (rolls) become stable at L < 0.45 (L > 0.45). At smaller ip we 
could compare this stability boundary between squares and rolls with results from Ref. 
and found excellent agreement. The oscillatory regime extends only up to L « 0.02. This, 



by the way, explains why oscillations could not be found in the experiments of Ref. |TTJ . For 
L > 0.02 squares transfer their stability directly to CRs when increasing r. 

Consider now the ip—r plane of Fig. |]b where stable 3D structures are predicted to appear 
at onset already for ip > 10~ 7 ||. For the smallest ip that we have analyzed, ip = 0.03, we 
found already the full sequence of squares, oscillations, and CRs with increasing r. With 
growing Soret coupling ip the stable existence range of each of these 3D structures that are 
peculiar to mixtures widens on cost of the 2D roll states. The onset of the latter is shifted 
upward to larger r: With increasing ip it requires higher thermal stress to generate the 
flow intensity for which advective mixing is strong enough to eliminate effectively the Soret 
induced concentration variations that ultimately cause and/or stabilize the 3D structures. 
Convection in well mixed binary fluids at large r is not much different from one-component 
fluids. Furthermore, our observation that the ratio of mean squared convective concentration 
and temperature variations that determine solutal and thermal buoyancy contributions, 
respectively, is large (small) compared to one at small (large) r sheds light on the oscillatory 
and/or stationary pattern competition at intermediate r: At small (large) r squares (rolls) 
are stably sustained by solutally (thermally) dominated buoyancy forces and in the crossover 
driving regime in between where neither squares nor rolls are stable there is oscillatory and/or 
stationary pattern competition. As an aside we mention that the advective concentration 
(heat) transport and the convective concentration (temperature) variation is slightly larger 
(smaller) in squares than in rolls. 

Studying the Prandtl number dependence in Fig. |c we found that oscillations disappear 
at cr < 0.2 and that rolls appear already slightly above r = 1. Increasing a the existence 
range of squares, oscillations, and CRs expands to higher r. We have also investigated 
parameters L, cr = (9(1) which apply to gas mixtures JTB|. Here typically rolls are stable 
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slightly above onset. Knowledge 17] of their stability domain is useful in relation to the 
competing spiral defect chaos that has recently been observed also in gas mixtures [16 . 
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FIGURES 



FIG. 1. Bifurcation diagram of x-roll intensity X and y-roll intensity Y versus r for squares 
(squares), rolls (circles), and crossrolls (triangles). Full (dashed) lines with filled (open) symbols 
denote stable (unstable) states. Full lines delimiting the hatched area mark maxima and minima of 
oscillations of X(t) and Y(t). A subharmonic bifurcation cascade occurs in the r-interval marked 
SC (width expanded for better visibility). For the parameters || tp = 0.23, L = 0.0045, a = 27 
used here stable squares and unstable rolls bifurcate with wave number k = ir out of the conductive 
state at r = 0.0124 while the critical Rayleigh number is r c = ^ = 0.00825. 

FIG. 2. Evolution of the oscillatory dynamics with increasing r from top to bottom. Parameters 
as in Fig. [I]. (a) and (b) are states located shortly above the Hopf threshold and right below the 
CR bifurcation, respectively, (c) - (f) are SC states. The right column contains phase space plots 
of X,Y versus X,Y for X(t),Y(t) shown in the left column with arbitrary common units versus 
reduced time n-. In the SC states n is the number of windings around each of the stable CR fixed 
points (triangles). Squares and circles indicate unstable square and roll states, respectively. 

FIG. 3. Grey-scale snapshots of topview shadowgraph intensity distribution I(x, y) evaluated 
with Eq. (||) for the oscillatory state of Fig.[2|b at times marked by A (y-rolls dominant), B (squares), 
and C ( x-rolls dominant). To elucidate the contributions from temperature and concentration we 
show in the bottom right corner the profiles of I T , I , and I = I T + I along the dotted line in B. 
Maximal up (down) flow is located in B at x = y ~ 0.5 (1.5). 

FIG. 4. Phase diagrams of stable convective states. The narrow SC domain is not shown. 
Parameters are ip = 0.23 (a, c), L = 0.0045 (b, c), and a = 27 (a,b). Fig. [l] shows for them a 
bifurcation diagram along the vertical dotted lines in (a, b). The ordinate scale is logarithmic for 
r < 1 to show the bifurcation threshold r s t a t{k = vr; tfi, L) of the conductive state (cond) which was 
determined with a shooting method. Phase boundaries in (a) between crossrolls and rolls from 
finite difference (full line) and Galerkin (dashed line) calculations are quantitatively unreliable 
where they disagree at small L combined with large r, i. e., when narrow concentration boundary 
layers occur in the fluid. 
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